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Abstract 

Background: Three pathogenicity islands, viz. SPI-1 (Salmonella pathogenicity island 1), SPI-2 (Salmonella 
pathogenicity island 2) and T6SS (Type VI Secretion System), present in the genome of Salmonella typhimurium 
have been implicated in the virulence of the pathogen. While the regulation of SPI-1 and SPI-2 (both encoding 
components of the Type III Secretion System - T3SS) are well understood, T6SS regulation is comparatively less 
studied. Interestingly, inter-connections among the regulatory elements of these three virulence determinants have 
also been suggested to be essential for successful infection. However, till date, an integrated view of gene 
regulation involving the regulators of these three secretion systems and their cross-talk is not available. 

Results: In the current study, relevant regulatory information available from literature have been integrated into a 
single Boolean network, which portrays the dynamics of T3SS (SPI-1 and SPI-2) and T6SS mediated virulence. Some 
additional regulatory interactions involving a two-component system response regulator YfhA have also been 
predicted and included in the Boolean network. These predictions are aimed at deciphering the effects of 
osmolarity on T6SS regulation, an aspect that has been suggested in earlier studies, but the mechanism of which 
was hitherto unknown. Simulation of the regulatory network was able to recreate in silico the experimentally 
observed sequential activation of SPI-1, SPI-2 and T6SS. 

Conclusions: The present study integrates relevant gene regulatory data (from literature and our prediction) into a 
single network, representing the cross-communication between T3SS (SPI-1 and SPI-2) and T6SS. This holistic view 
of regulatory interactions is expected to improve the current understanding of pathogenesis of 5. typhimurium. 

Keywords: Salmonella typhimurium, Salmonella pathogenicity island 1 (SPI-1), SPI-2, Type VI Secretion System (T6SS), 
Boolean modeling, Cross-talk network 



Background 

Salmonella enterica serovar typhimurium is one of the 
causative agents of gastroenteritis in humans. The infec- 
tion caused by this pathogen, known as salmonellosis, 
leads to severe fever, diarrhea and even death in 
immune-compromised patients. Such infections, al- 
though more prevalent in developing countries (with 
most reported cases in Asia), continue to be a global 
threat with a significant number of non-typhoid salmon- 
ellosis incidences reported from industrialized countries. 
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Unraveling the mechanism of infection by the pathogen 
has been a major focus for researchers working in the 
field of enteric diseases. Two Type III Secretion Systems 
(T3SS) have been primarily implicated to be associated 
with the virulence of S. typhimurium. The genes encod- 
ing the T3SS components are located in two pathogen- 
icity islands, namely, SPI-1 and SPI-2 (Salmonella 
pathogenicity islands 1 and 2), present in the genome of 
S. typhimurium [1,2]. SPI-1 and SPI-2 differ significantly 
in terms of the types of genes present in them, their 
regulation, and the environment in which they are 
expressed. The initial stages of infection by the pathogen 
are controlled by the genes of SPI-1, which enable the 
pathogen to invade the epithelial cells of the intestine 
and induce an inflammatory response [1]. On the other 
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hand, once the pathogen invades the host macrophages 
and proliferates intra- cellularly, genes corresponding to 
SPI-2 get expressed [3,4]. Interaction among the major 
regulators of these two secretion systems (encoded by 
SPI-1 and SPI-2) have also been shown to play crucial 
roles in mediating virulence [5,6]. 

In addition to T3SS, the Type VI Secretion System 
(T6SS) has also been reported to be functionally active in 
S. typhimurium [7]. T6SS has been suggested to play a sig- 
nificant role in later phases of infection. Although the 
genes constituting T6SS do not participate directly in 
pathogenesis, they apparently enhance the efficiency of 
colonization/ infection [8]. It has been proposed that, T6SS 
genes in S. typhimurium get expressed in the host macro- 
phage at a late stage of infection when the levels of SPI-2 
proteins are reduced [9]. Intriguingly, this inverse regula- 
tion between SPI-2 and T6SS has been implicated in help- 
ing the pathogens persistence inside macrophage by 
preventing its over-growth and subsequent rapid death of 
the host cell [9]. In summary, a wide range of experimen- 
tal evidences suggest that a sequential activation of SPI-1, 
SPI-2 and T6SS are required during the progression of the 
chronic infection. A schematic representation of different 
stages of infection mediated by these secretion systems of 
the pathogen is depicted in Figure 1. The different envi- 
ronments experienced by the pathogen in course of infec- 
tion have also been included in the figure. 

Previous studies on secretion systems in S. typhimurium 
indicate existence of a tightly controlled gene regulatory 



network mediating cross-talk among SPI-1, SPI-2 and 
T6SS [1,4,9]. In order to achieve an opportune and com- 
plex regulatory control over the target gene expression, 
gene regulatory networks are expected to rely on sensing 
environmental cues using sensor- regulator systems as well 
as on interactions among different transcription regula- 
tors. Although, these two aspects have been widely studied 
for the network representing interconnections between 
regulators of SPI-1 and SPI-2 in Salmonella, the cross-talk 
network between T3SS and T6SS regulators is not yet well 
characterized. For example, several small-scale network 
models have been proposed for the gene regulation of the 
two T3SSs as well as their interconnections [1,5,10-14] 
mediated through transcription regulators (like HilD) and 
two-component systems (like PhoP-PhoQ). On the other 
hand, relatively fewer interactions related to the communi- 
cation between T3SS and T6SS have been reported. For 
instance, interaction between two-component systems 
PhoP-PhoQ (involved in regulation of both SPI-1 and SPI- 
2) and PmrA-PmrB (involved in regulation of T6SS) have 
been reported by Soncini and Groisman [15]. An indirect 
regulation of RcsB (another major regulator of T6SS) by 
PhoP, mediated via the response regulator MviA, has also 
been identified in a subsequent study [16]. One of the 
inner membrane components of T6SS, namely SciS, has 
also been shown to be negatively regulated by SsrB (tran- 
scription regulator involved in SPI-2 activation) [9]. How- 
ever these interactions, when viewed in isolation, cannot 
explain the sequential activation of the secretion systems. 
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Figure 1 Salmonella infection. Sequential activation of different secretion systems of Salmonella during infection. The levels of key 
environmental factors sensed by the pathogen - first in the distal gut lumen, and subsequently in the intra-macrophageal environment, 
are mentioned. 
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Unlike for T3SSs, the sensing of environmental signals 
driving T6SS regulation is not well-defined. For example, 
although osmolarity has been suspected to have effects on 
the transcription regulator RcsB [17,18], the exact under- 
lying mechanism has not been elucidated. Therefore, iden- 
tification as well as integration of the missing interactions 
and environmental triggers would help in better represen- 
tation of the cross-talk between T3SS and T6SS, there- 
by strengthening our understanding of S. typhimurium 
pathogenicity. 

The current study focuses on integrating relevant 
regulatory information associated with SPI-1, SPI-2 and 
T6SS into a single network and subsequently designing 
simple mathematical model(s) to gain a systems-level 
view of gene regulation driving the virulence of S. 
typhy murium, 'Boolean modeling', particularly useful in 
cases where there is a scarcity of experimental data [18], 
has been used in the present study. In Boolean models, 
each node in a network can attain a binary state, which 
is either 'one' (ON) or zero' (OFF), indicating the pres- 
ence/activation or absence/inactivation of the element 
represented by the node. At any given time point, the 
state of each node is determined by the states of its 
neighbors (interacting with it) at the previous time point. 
The rule governing the change of state of any node is 
thus framed depending on the nature of the interactions. 
In other words, the rules can be logical statements in- 
volving the states of corresponding regulatory nodes 
which are joined using logical operators like 'AND', 'OR', 
'NOT; etc. In spite of the simplicity of its logical frame- 
work, Boolean models have been successfully employed 
to study and predict the dynamic behavior of gene regu- 
latory networks [19-25]. Furthermore, the sufficiency of 
the interactions in a gene regulatory network can be 
verified through Boolean modeling, based on which novel 
interactions may be proposed to fill potential gaps [19]. In 
view of this, a Boolean model of the gene regulatory net- 
work of SPI-1, SPI-2 and T6SS of S. typhymurium was 
constructed. An attempt was made to identify previously 
unreported response regulators (or transcription factors) 
that may be participating in activation of T6SS. Subse- 
quently, simulation(s) were run with the model in order to 
verify whether the interactions integrated in this network 
(including the ones predicted in this study) could mimic 
the experimentally observed sequential activation of the 
three secretion systems. 

Results and discussion 

Compilation of relevant regulatory interactions into a 
Boolean model 

The first step towards modeling a gene regulatory net- 
work involves compilation of available interaction data 
from literature. An exhaustive literature survey was 
performed to identify relevant genes and significant 



regulatory interactions that are involved in activation of 
SPI-1, SPI-2 and T6SS of Salmonella enterica serovar 
typhimurium. Since environmental factors (such as osmo- 
larity, metal ion concentrations, availability of sugars, etc.), 
encountered by the pathogen inside the host during differ- 
ent stages of infection, affect the regulation of these genes, 
available information about them and their regulatory 
roles were also compiled. These data provided the founda- 
tion for construction of the Boolean model representing 
the regulation and cross-talk network of Type III and 
Type VI Secretion Systems in S. typhimurium. 

For constructing the Boolean model, the gene regulatory 
interactions were represented as a 'directed graph'. A sche- 
matic diagram of the network is shown in Figure 2. Each 
node in this graph symbolized either one gene or an en- 
semble of genes working in tandem (e.g. two-component 
systems, components of SPI-1, SPI-2 and T6SS). Various 
environmental factors were also represented by nodes. 
Each (directed) edge in this graph depicted either a posi- 
tive or a negative regulatory effect exerted by a particular 
node on its neighbor or on itself (in case of self-regulation). 
The genes included in the network are listed in Additional 
file 1. The gene regulatory cascade leading to virulence is 
expected to be triggered on receiving certain environmental 
cues or at certain stages of growth/development of the 
pathogen. Therefore, except such environmental factors, no 
orphan nodes' (nodes which lack incoming regulatory con- 
nections) were kept in the network. Thus, the network was 
constructed to closely resemble the in vivo regulatory inter- 
actions involving the virulence factors of S. typhimurium 
and various environmental factors at the sites of infection 
inside human body. 

Transcription factor binding sites in the neighborhood of 
T6SS genes 

Since the regulation of T6SS in Salmonella is not well 
studied, it may be anticipated that a regulatory network 
constructed with the existing experimental data would 
lack the expected repertoire of activating factors as well 
as environmental sensors that are required to trigger 
T6SS. For example, although previous studies have sug- 
gested the role of osmolarity in regulation of one of the 
principal activators of T6SS, namely RcsB [17,18], the 
exact mechanism, the sensor-regulator system as well as 
other intermediates involved in this activation cascade, 
have not been identified yet in S. typhimurium. It was 
therefore alluring to investigate whether the two regula- 
tors of T6SS (viz., RcsB and PmrA) included in the com- 
piled network of interactions can be connected to any 
additional regulatory cascades sensing the environment, 
especially osmolarity. The mechanism of osmoregulation 
of SPI-2 by the two-component system OmpR-EnvZ is 
well established [26]. Intriguingly, a study on two- 
component systems in Eschericia coli has shown that 
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Figure 2 SPI-1, SPI-2, and T6SS regulatory network. (A) Constructed gene regulatory network controlling secretion systems of S. typhimurium. 
A Boolean model of this regulatory network was designed by defining appropriate logical rules for the depicted interactions. (B) The initial states 
defined for simulating the Boolean model under two different environments, viz., distal intestinal lumen and the macrophage. 



several sensor kinases have the ability to pair up with 
multiple non-cognate response regulators [27], thereby 
facilitating cross -talk between various sensor-regulator 
pathways. It has also been observed that the osmolarity 
sensor EnvZ, whose cognate response regulator is 
OmpR, can additionally phosphorylate a non-cognate re- 
sponse regulator YfhA in E. coli. YhfA is normally known 
to be a cognate response regulator of the YfhK-YfhA 
two-component system [28], which gets activated during 
stationary phase of growth (a condition associated with 
T6SS activation in S. typhimurium). BLAST similarity 
searches indicated that the protein STM2562 of S. 
typhimurium is a close homolog of the E.coli YfhA (95% 
identity, 100% coverage). STM2562 has also been anno- 
tated to have a similar function in the NCBI database. 
Considering that the osmolarity sensor EnvZ is not 
known to have any other non-cognate response regula- 
tor, YfhA was further investigated for its possible role in 
osmoregulation of T6SS. Since most response regulators 
control gene expression at the transcription level by 



binding to the DNA, the neighborhood of T6SS genes 
were searched for YfhA binding sites. In addition, 
the Tfsitescan server (http://www.ifti.org/Tfsitescan/) was 
employed to identify the presence of any other transcrip- 
tion factor binding motifs near the T6SS genes. 

(A) YfhA binding sites 

A search for YfhA binding sites with a previously reported 
18 base pair (bp) binding motif [28] revealed three prob- 
able hits in the upstream of the rcsB gene, when a mis- 
match of 1 bp was allowed (Additional file 2). Allowing a 
mismatch upto 2 bp led to the identification of four prob- 
able binding sites within the ORF of rcsB in addition to 
multiple hits in the upstream region (Additional file 2). 

(B) FlrC binding sites 

Multiple probable transcription factor binding mo- 
tifs were detected using the 'Tfsitescan server (see 
'Methods') in the upstream of T6SS components (SciS, 
VrgS and SciG) and their regulators (RscB and PmrA) in 
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S. typhimurium. These are listed in Additional file 3. 
One of the detected binding sites corresponded to the 
transcription factor FlrC of Vibrio cholerae. The up- 
stream regions of rcsB and sciS were found to contain 
copies of a six nucleotide long motif (CGGCAA) identi- 
cal to the FlrC binding site [29]. Further searches with 
Tfsitescan resulted in detection of additional FlrC bind- 
ing sites in the ORFs of rcsB and sciS. Details of the FlrC 
binding sites found within the ORFs and upstream re- 
gions of sciS and rcsB are listed in Additional file 4. FlrC 
has been reported to function in flagellar biogenesis and 
intestinal colonization of V. cholerae [29]. Moreover, the 
flagellar system has been shown to play a significant role 
in regulation of T6SS in V. cholerae and have also been 
implicated in the regulation of SPI-1 in S. typhimurium 
in earlier studies [30,31]. A homology search in S. 
typhimurium indicated the protein YfhA to be the clos- 
est homolog of FlrC (of V. cholerae) with a sequence 
similarity of 62%, strengthening the possibility of YfhAs 
involvement in regulation of T6SS in S. typhimurium. 

It may be noted here that Tfsitescan could not identify 
any binding sites for YfhA, probably due to absence of 
this information in its backend database. Also, only a 
few of the other binding sites identified could be associ- 
ated with transcription factors directly involved in viru- 
lence related pathways (Additional file 3). 

Structural homology between YfhA of S. typhimurium and 
FlrC of V. cholerae 

A structural comparison of Vibrio FlrC and Salmonella 
YfhA (having 62% sequence similarity) was performed to 
check the extent of their structural homology, which 
would imply a similar function. Since no crystallographic 
structures of either of the proteins were available, the 
comparison was done by building homology models of 
these two proteins. One of the transcription regulator 
proteins belonging to the NtrC family from Aquifex 
aeolicus (PDB ID: lny5B), having the best matches with 
each of these two proteins, was considered as the tem- 
plate for building their models. When the modeled 
structures of FlrC and YfhA were superposed, the RMSD 
was observed to be 1.32 A (Additional file 5). A 
segment-matching approach [32] was subsequently 
adopted for further evaluating the anticipated structural 
similarity. FlrC and YfhA are reported as response regu- 
lators of two-component systems, and are known to par- 
ticipate in sigma-54 dependent transcription [28,29]. 
Such response regulators are known to have three dis- 
tinct functional domains, viz., a response regulator re- 
ceiver (REC) domain, a sigma-factor interaction (AAA) 
domain, and a helix-turn-helix DNA binding (HTH) do- 
main. The sequence segments corresponding to these 
three distinct domains of FlrC and YfhA were identified 
by searching against the Pfam database [33] (Additional 



file 6). These domains/segments were then modeled sep- 
arately using the best available templates (in PDB) corre- 
sponding to each of these domains (Additional file 7). 
Subsequently, pair-wise structural superposition was 
performed using the models for each of the domains of 
FlrC and YfhA (Additional file 7). Figure 3 shows the in- 
dividual and superposed structures of these domains 
from both the proteins. The RMSD between the respect- 
ive domains were observed to be 0.09 A for the REC do- 
mains, 0.79 A for the HTH domains and 1.47 A for the 
AAA domains. These results suggest that the YfhA pro- 
tein in S. typhimurium shares significant similarities in 
structure with the FlrC protein in V. cholerae, and may 
therefore participate in similar functional activities. 

Probable role of YfhA in osmoregulation of T6SS 

The presence of several binding sites of YfhA (and its 
homolog FlrC) in the upstream regions and ORFs of rcsB 
and sciS, coupled to the fact that YfhA is also a non- 
cognate response regulator of EnvZ, strongly indicates the 
involvement of YfhA in osmoregulation of T6SS. This 
mechanism is probably similar to that of EnvZ-OmpR me- 
diated osmoregulation of the SsrA-SsrB system [34]. It has 
been reported that, high osmolarity results in increased 
phosphorylation of OmpR, a cognate response regulator 
of EnvZ [34]. In low osmolarity conditions, OmpR and/or 
OmpR-P (phosphorylayted OmpR present in low levels) 
can bind to the upstream regions of ssrA and ssrB genes. 
OmpR-P, however, have a higher DNA binding affinity as 
compared to that of OmpR, and when in abundance (as in 
high osmolarity conditions) can bind to the ORFs of ssrA 
and ssrB, thereby blocking transcription of these two genes 
[34]. In the present study, rcsB and sciS were also pre- 
dicted to have binding sites of YfhA (and its homolog 
FlrC) within their ORFs (Additional file 4). Hence, it is 
likely that there exists a similar mechanism of phosphory- 
lated YfhA (YfhA-P) mediated regulation of sciS and rcsB. 
In view of this, it may be hypothesized that, under low 
levels of osmolarity (typical of an intra-macrophageal en- 
vironment), unphosphorylated YfhA acts as a sigma-54 
dependent transcription factor for the genes sciS and rcsB, 
and is responsible for activating the T6SS expression. On 
the other hand, once the response regulator YfhA gets 
phosphorylated by EnvZ under high levels of osmolarity, it 
is likely that it additionally binds to the ORFs of sciS and 
rcsB, thus limiting the transcription of these two genes. 
The missing links connecting osmolarity and RcsB in the 
Boolean model were imputed based on these assumptions. 

Sequential activation of SPI-1, SPI-2 and T6SS through 
simulation of Boolean models 

Simulation of the Boolean model was performed in 
order to study the dynamic behavior of the regulatory 
network and to verify the fitness of the predicted 
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Figure 3 Structural homology of FIrC and YfhA. Model structures of REC, AAA, and HTH domains of FIrC and YfhA proteins. The bottom panel 
shows pair-wise superposition of corresponding domains and their RMSD. 



interactions involving YfhA in the model The Boolean 
model was simulated under both synchronous' and 
asynchronous' modes using "GINsim" (http://ginsim. 
org/home) [35]. The results of the synchronous simula- 
tions (state transitions for every node in the network) 
are provided in Table 1. While simulating the model 
(allowing synchronous updates) under intestinal envir- 
onment, the model traversed through seven states before 
reaching a steady state with SPI-1 turned 'ON' and other 
two secretion systems (SPI-2 and T6SS) turned 'OFF'. 
Subsequently, when the model was simulated under 
intra-macrophageal environment, it converged to a 
steady state having T6SS turned 'ON' and the two T3SSs 
turned 'OFF', after traversing five intermediate steps. 
However, during the second simulation, SPI-2 was ob- 
served to get activated and then repressed in two inter- 
mediate steps before arriving at the steady state. The 
results of both the simulations taken together represent 
a sequential activation of SPI-1, SPI-2 and T6SS, which 
corroborates well with those observed experimentally 
[5,9]. 

Synchronous simulations, however have a limitation, 
considering that for the same initial condition, it would al- 
ways progress through the same path (transition states) 
leading to the same stable state/cycle. To investigate other 
possible stable states (and intermediate steps), multiple 
synchronous simulations were performed wherein only 



the nodes corresponding to the environmental factors 
were initialized (fixed) as in earlier simulations. The regu- 
latory nodes for these simulations were set to random ini- 
tial states (in all possible combinations). Furthermore, 
simulations were also performed using an asynchronous 
update scheme. Adopting an asynchronous update scheme 
makes the outcomes of the simulation stochastic and al- 
lows to find out multiple steady states, in case they exist. 
All these simulations led to the same two steady states, 
which were identical to the steady states obtained through 
the synchronous simulations performed earlier (Table 1). 
The transition probabilities between different states of the 
system, during asynchronous simulations, were also 
analysed for the intra-macrophageal environment. The re- 
sults from the asynchronous simulations (Additional 
file 8) indicate that the system has a high chance of 
progressing through paths which correspond to the pro- 
posed sequential activation. These results suggest that the 
Boolean model constructed in the present study fairly cap- 
tures the in vivo regulatory mechanism controlling the 
cross-talk between T3SS and T6SS in S. typhimurium. 

An intriguing behavior of the system under study corre- 
sponds to the activation of SPI-2 and T6SS under similar 
environmental conditions inside the macrophages at dif- 
ferent stages of the infection cycle. Although, previous 
studies have suggested a relation between low concentra- 
tion of SPI-2-secreted proteins and activation of T6SS [9], 
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(B) 



Nodes 



State Glucose Iron Magnesium Calcium Osmolarity Stationary MIc HilE HilD HilC RtsA HilA Ihf SirA- CsrBC CsrA Fur H-NS PhoP SlyA SsrAB Fis EnvZ OmpR YfhA MviA RcsB PmrA SciS VrgS SciG SPI -1 SPI -2 T6SS 
Transition phase BarA 
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The rows represent the dynamic behavior of the network (constituting 34 nodes) under study. '0' represents inactive state and T represents active state. Row 1 represents the values specified as the initial states for 
simulation. A node's transit from an inactive state to an active state (0 — > 1) or vice-versa (1 — > 0) has been marked with asterisks (*). (A) Simulation considering the distal intestinal lumen environment, which is the 
site of infection by 5. typhimurium. The 8th row corresponds to the steady state attained by the network (representing an activated SPI-1) (B) Simulation considering the intra-macrophageal environment faced by S. 
typhimurium after it has been engulfed by the macrophage. The 6th row corresponds to the steady state attained by the network (representing an activated T6SS). 
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the mechanism of the regulation is still unclear. From the 
network constructed in the current study, it can be in- 
ferred that several antagonistic connections, such as 
SsrAB to SciS and MviA to RcsB, may be responsible for 
keeping T6SS in an inactive state while SPI-2 remains ac- 
tive. Analysis of the state transition table obtained through 
model simulation shows that, a major inner membrane 
component of T6SS (SciS) becomes active only when 
SsrAB is turned off, thus making SPI-2 inactive. SsrAB 
and SciS are otherwise observed to sense similar environ- 
mental signals (like low osmolarity, low levels of magne- 
sium and calcium) through different routes (direct or 
indirect), which allows activation of both SPI-2 and T6SS 
in an intra-macrophageal environment. These results sug- 
gest that, coordinated activation of SsrAB and SciS, prob- 
ably plays a crucial role in sequential activation of SPI-2 
and T6SS respectively. 

To further understand the role of YfhA in the model, a 
simulation was run under synchronous mode with a (sim- 
ulated) mutant of this protein, keeping the initial condi- 
tions unaltered. The mutation was introduced using the 
mutant simulation' feature provided in 'GINsiml which re- 
stricted the YfhA node to an inactive (zero) state through- 
out the entire simulation. Although including this mutant' 
in the simulation didn't affect the sequential activation 
and repression of SPI-1 and SPI-2, the model failed to 
reach a state with T6SS turned 'ON' (Additional file 9), in- 
dicating the importance of YfhA in the network. It may 
however be noted that this outcome also relies on the 
underlying assumption (while defining the logical rules) 
that sensing osmolarity and magnesium concentrations 
are equally essential for an opportune T6SS activation. 

Conclusions 

A Boolean modeling approach has been utilized in the 
present study to understand the dynamics of the secre- 
tion systems mediating virulence in S. typhimurium. In 
order to unravel the underlying molecular mechanism of 
virulence of this pathogen, two Salmonella pathogenicity 
islands (SPI-1 and SPI-2) have been the focus of many 
earlier studies. Existing literature provides a rich reper- 
toire of gene regulatory information and indicates in- 
volvement of a complex cross-talk among various 
regulatory elements of SPI-1, SPI-2 and T6SS in the 
pathogen's virulence. In order to obtain a holistic view of 
the regulatory interplay, the present study integrates 
relevant gene regulatory data into a single network, 
representing the cross-communication between T3SS 
(SPI-1 and SPI-2) and T6SS. This comprehensive view is 
expected to improve the current understanding of 
pathogenesis of S. typhimurium. In addition, involve- 
ment of new components and interactions has been pre- 
dicted in an effort to fill in some gaps in our current 
knowledge of this network. A major challenge towards 



modeling the cross-talk network under a Boolean frame- 
work pertains to properly capturing the intricacies of 
gene regulation events. This relates to accurately defin- 
ing the logical rules to ensure the correctness of the 
model, and to use only the environmental factors (that 
are critical for triggering any gene regulatory cascade) as 
input nodes. The logical rules have been defined to max- 
imally reflect the experimentally observed modes of gene 
regulation involved in the network under study (de- 
scribed in the 'Methods' section). The environmental 
factors have been connected to the network in order to 
mimic different extra- and intra-cellular environments 
experienced by the pathogen during various stages of in- 
fection. Boolean modeling of the cross-talk network and 
subsequent simulation could reproduce the experimen- 
tally observed sequential activation of SPI-1, SPI-2 and 
T6SS, indicating that, the model could efficiently capture 
the underlying regulatory mechanism driving the pro- 
gression of infection in vivo. 

An interesting outcome of the current study is the 
identification of YfhA as a potential transcriptional regu- 
lator involved in the sub-network corresponding to 
T6SS. The constructed sub-network for T6SS regulation, 
based on the information from literature, contains only 
two regulatory elements (RcsB and PmrA), which have 
been reported as the major regulators of T6SS genes 
[36]. Additionally, the sub-network could be connected 
to only one of the environmental factors (namely, mag- 
nesium) through the transcriptional regulator PhoP [37]. 
Although the suggested role of osmolarity in regulation 
of RcsB has been reported earlier [17,18], the sensor- 
regulator cascade remains to be identified experimen- 
tally. Imputing YfhA as a missing link would therefore 
enrich the model, considering that sensing both the en- 
vironmental factors (low osmolarity, low magnesium) is 
likely to result in a stronger signal for the T6SS activa- 
tion cascade, than sensing either one of them. However, 
such an effect would be difficult to capture through a 
model based on Boolean-logic, and could probably be bet- 
ter investigated with an ODE-model. Introduction of YfhA 
in the T6SS sub-network also results in formation of a 
feed-forward loop involving YfhA, RcsB and SciS. A recent 
investigation by Le and Kwon (2013) has indicated that 
such coherent feed-forward loops (FFLs) improve network 
robustness against update-rule perturbations [38]. Al- 
though, validating the effect of this FFL in improving net- 
work robustness is outside the scope of this work, the 
inclusion of YfhA (and the FFL) probably makes the T6SS 
network better equipped to sense and respond to the en- 
vironmental cues. 

In the current study potential binding sites of YfhA, a 
two-component system response regulator, has been 
identified in the upstream region and ORF of the rscB 
gene (one of the major regulators of T6SS). YfhA has 
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also been shown to act as a non-cognate response regu- 
lator of the histidine kinase EnvZ [27]. Since EnvZ is a 
known osmolarity sensor for SPI-2 genes [26], it is likely 
that YfhA regulates T6SS activation (through RcsB) de- 
pending on the osmolarity sensed by EnvZ. Binding mo- 
tifs recognized by FlrC, a transcriptional regulator found 
in V. cholerae, has also been detected in the upstream 
regions and within the ORFs of rcsB and sciS (an inner 
membrane component of T6SS). Cooperative binding of 
FlrC near the downstream of transcription start sites of 
two flagellar genes (flaA and flgK) has been shown to en- 
hance intestinal colonization of V. cholerae [29]. Though 
FlrC is not found in S. typhimurium, analyses performed 
on the protein sequences as well as their structures sug- 
gest that the protein YfhA (from S. typhimurium) is a 
close homolog of FlrC. It may be noted that, the binding 
sites predicted in this study are based on the motifs recog- 
nized by two YfhA homologs from E. coli and V. cholerae. 
The motifs recognized by the YfhA of S. typhimurium may 
be different from the ones screened for in this study, 
which can only be confirmed with further experimental 
validation. However, these predictions coupled with the in- 
formation available from literature suggest that inclusion 
of YfhA in the model is able to bridge the gap pertaining 
to the anticipated, but hitherto unknown, osmoregulatory 
cascade controlling T6SS activation. Furthermore, simula- 
tion of the Boolean model with a mutated' YfhA did not 
result in T6SS activation, indicating the probable crucial 
role of YfhA in the constructed network. 

The current study presents a comprehensive model 
representing cross-talk among the regulatory elements of 
the three major pathogenicity islands in S. typhimurum 
and predicts probable additional regulatory links to impute 
certain gaps in the network. Mapping of the dynamic be- 
havior of the cross-talk network suggests the role of coor- 
dinated regulation of SsrAB and SciS in successive 
activation of SPI-2 and T6SS, rather than simultaneous ac- 
tivation. Results obtained in this study are expected to 
benefit researchers in understanding the sequential activa- 
tion of secretion systems in S. typhimurium and the pro- 
gression of infection. The integrated view of the gene 
regulatory network presented here will aid in future studies 
on S. typhimurium infection. Further experimental valid- 
ation of these findings would strengthen our understanding 
of S. typhimurium pathogenesis. In addition, future efforts 
towards understanding how host proteins interact with the 
secreted virulence factors over the course of infection is 
likely to provide a broader overview. 

Methods 

Search for novel regulatory cascades for T6SS activation 
(A) Search for YfhA binding sites 

A previous study on E. coli had indicated a 18 bp motif 
recognized by YfhA [28]. Multiple copies of this motif 



[TGTCN(10)GACA] has been shown to be present up- 
stream of the glmY gene in E. coli which is regulated by 
YfhA. Presence of this motif was searched in the upstream 
regions of the major components of T6SS, namely, sciS, 
vrgS, sciG, as well as the direct regulators of T6SS compo- 
nents, namely, rcsB and pmrA. In cases where a gene was 
preceded by other genes in its respective transcription unit 
(operon), the search was performed in the upstream re- 
gions of the corresponding operon. Information about pu- 
tative transcription units were retrieved from DOOR 
(Database of prokaryotic operons) [39]. 

(B) Search for other transcription factor binding sites 

The upstream regions of the T6SS components and their 
direct regulators were also screened for the presence of 
other Transcription Factor (TF) binding sites. The tool 
"Tfsitescan" (http://www.ifti.org/Tfsitescan/) [40] was used 
to detect such TF binding sites. Subsequently, the list of 
transcription factors reported by "Tfsitescan" were filtered 
based on their functions (such as their association to viru- 
lence) in order to narrow down on probable regulatory 
candidates which could also have connections to the antic- 
ipated EnvZ-YfhA mediated osmolarity sensing. 

(C) Structural study of YfhA and FlrC 

The search for TF-binding sites indicated several occur- 
rences of the TF-binding motif corresponding to one of 
the flagellar regulatory proteins (FlrC) of V. cholerae (see 
'Results and Discussion'). Sequence similarity searches 
also asserted FlrC to be the closest homolog of the re- 
sponse regulator YfhA from S. typhimurium. The extent 
of structural relatedness (implying a functional similar- 
ity) between these homologs was compared using a 
homology modeling approach. The model structures 
were generated with the help of Swiss Model server [41] 
by providing appropriate template structures (see 'Re- 
sults and Discussion'). For checking the structural simi- 
larity, the homology models were superposed with the 
help of the freely available software 'FAST Alignment 
and Search Tool' [42] and RMSD value (s) were noted. 

Boolean model construction 

The network of interactions was modeled under Boolean/ 
logical framework using the tool 'GINsim' (http://ginsim. 
org/home), which provides a user-friendly interface for 
simulation and efficient analysis of regulatory graphs [35] . 
Model building using GINsim involved the following three 
steps, namely, (i) construction of a regulatory network 
consisting of nodes and edges, (ii) defining logical rules for 
each node, and (iii) defining basal and maximum levels of 
expression for each node. The network was constructed 
using the experimentally verified interactions obtained 
from literature. In addition, the predicted regulatory cas- 
cade involving YfhA (see 'Results and Discussion'), was 
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also included. The network contained 34 nodes 
representing genes, their regulators, as well as environ- 
mental factors, and 65 edges representing positive or 
negative regulations. The logical rules for each node were 
defined based on the modes of regulation of the node by 
one or more regulators. The logical rules and correspond- 
ing references to literature are provided in Additional 
file 10. In case of a node having multiple regulators (positive 
and/or negative), AND' or 'OR' operators were used in the 
rule based on the relative effects of the regulators on the 
activation/repression of the particular node. For example, 
HilA, the main regulator directly activating SPI-1 genes, 
has been reported to get sufficiently expressed in the pres- 
ence of HilD along with HilC or RtsA [11]. Hence the rule 
defined for HilA was "HilD AND (HilC OR RtsA)". For 
cases, where adequate information regarding the relative 
effects of the regulators were not available, the 'OR' logic 
was applied. One such case was observed for HilD, a 
major transcriptional regulator involved in expression of 
SPI-1 genes. In this case, HilD was subjected to positive 
regulation by four nodes (SirA-BirA, HilC, RtsA and Fur) 
along with self activation [11,14,43,44]. Negative regula- 
tions were included in the rules (wherever reported) with 
an AND NOT operator, since all negative regulations 
were reported to be significantly repressive (except in case 
of H-NS). The rule defined for H-NS, a gene involved in 
nucleoid organization, differed from this standard assump- 
tion. H-NS was observed to have five negative regulators 
(HilD, PhoP, SlyA, Fur and IHF), either direct or indirect, 
relative effects of which were not available from literature. 
While H-NS has been reported to act as a repressor for a 
wide range of genes, H-NS mediated repressions have 
been shown to be counter-acted by different regulators. 
Such examples include counter-effect of HilD on H-NS 
mediated silencing of ssrA and ssrB [5] and counter-effect 
of IHF on H-NS mediated repression of hilA [45]. To ac- 
count for such multi-factorial effects, the rule for H-NS 
was defined such that H-NS was turned 'OFF' only in the 
presence of all five repressors (HilD, PhoP, SlyA, IHF and 
Fur). In other cases wherein one or more of the repressors 
of H-NS were absent, the rules were appropriately defined 
to consider the counter-effects on H-NS mediated repres- 
sion (without turning 'OFF the H-NS node itself). 

The basal and maximum levels for the nodes corre- 
sponding to the environmental factors were defined based 
on the information obtained from the existing literature. 
The reported availability of these factors in the distal intes- 
tinal lumen (site of initial stage of infection) and in macro- 
phages (site of late stages of infection) infected with S. 
typhimurium were considered for the model. Glucose has 
been shown to mostly get absorbed in the small intestine, 
limiting its availability in the distal intestine [46]. On the 
other hand, inside the macrophage, the pathogen has been 
reported to require glucose for causing successful 



infection [47]. Osmolarity, magnesium and calcium have 
been found to be significantly higher in distal intestine as 
compared to those in the macrophage [14,48-50]. Thus to 
accommodate for the changes in environment, the basal 
and maximum levels were kept as zero' and one' respect- 
ively for 'Glucose! 'Osmolarity 'Magnesium' and 'Calcium'. 
The concentration of iron has been shown to remain com- 
paratively high in the intestine due to absorption of dietary 
iron [14]. Also during later stages of infection, high iron 
concentration has been reported to aid in bacterial sur- 
vival and multiplication inside the macrophage [37]. Thus, 
to account for high iron concentration both during initial 
(inside distal intestine) and later phases (inside macro- 
phage) of infection, the basal and maximum values for 
'Iron' were defined as 'one'. Similar values were assigned 
for the node 'Stationary phase! as, under oxygen-limiting 
condition (the site of pathogenic colonization inside hu- 
man body), S. typhimurium in stationary phase has been 
shown to exhibit optimal invasiveness [51,52]. Among 
other nodes, the basal and maximum levels of the nodes 
having at least one positive incoming connection (e.g., 
HilD, Fur, RcsB, etc.) were defined as 'zero' and 'one! re- 
spectively. For the nodes having only negative incoming 
connections (e.g., HilE, CsrA, H-NS, etc.), both the levels 
were made as 'one' assuming that, they are constitutively 
expressed in absence of any repressors. The Boolean 
model was thus constructed by integrating all the relevant 
experimental observations. It may however be noted that 
the environmental factors considered in the current study 
only provides a simplistic view of the in-vivo conditions en- 
countered by the pathogen. There can be several other fac- 
tors associated with the virulence mechanism. For 
example, phagocytic free radicals (like hydroxyl radical 
HO-) have been reported to cause DNA damage. In re- 
sponse to such DNA damage, intra-cellular Salmonella has 
been reported to induce SOS response, which mainly in- 
cludes differential expression of genes like recA, uvrABY, 
lexA, sulA, sodA, sodB, etc. [53,54]. Previous studies have 
also suggested an SPI-2 dependent mechanism in S. 
typhimurium which inhibits the production of reactive 
oxygen species [55]. However, no reported regulatory inter- 
action^) could be found (from existing literature) which 
indicate any direct role of SOS response (or other such fac- 
tors) in activating any of the three secretion systems con- 
sidered in this study. As more experimental data gets 
available in future, it should be possible to incorporate sev- 
eral such environmental factors and response mechanisms 
to extend the proposed model. 

Simulation of the Boolean model 

The Boolean model was simulated with GINsim [35], in 
order to draw inferences about the dynamic behavior as 
well as the steady states of the system. Simulations were 
performed using both 'synchronous' and 'asynchronous' 
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modes. 'Synchronous' mode offers a deterministic up- 
date schema, wherein each node is updated once in a 
single time step, leading to change in the network state 
at the end of each update cycle [56]. Adopting an asyn- 
chronous' update scheme allows updation of any ran- 
domly selected node at a time step, thereby making the 
outcomes stochastic. While performing the simulation 
runs, only the environmental factors were considered as 
input nodes. The environmental factors included in the 
current study were osmolarity, stationary phase of 
growth as well as availability of glucose, magnesium, cal- 
cium and iron. 

Activation of SPI-1 is known to lead to the invasion of in- 
testinal epithelial cells by S. typhimurium, and subsequent 
engulfment of the pathogen by macrophages. Thus, in 
order to reflect this change in environment experienced by 
the pathogen, the model under study was simulated with 
two different sets of initial conditions. These conditions 
corresponded to the environments in the distal intestinal 
lumen and in the macrophage. While simulating the model 
for the distal intestinal lumen environment, the initial state 
for 'Glucose' was set at zero! and other environmental fac- 
tors were set at one! according to their expected availability 
at the beginning of the infection cycle (described in the pre- 
vious section). For all other nodes, various possible combi- 
nations of zero' and one' were considered as initial states 
for different rounds of simulation, to witness the behavior 
of the system constrained only by the environmental cues 
and the defined logical rules. Similarly, while simulating the 
model for intra-macrophageal conditions, the values for 
the environmental factors 'Glucose,' Iron' and 'Stationary 
Phase' were set as 'one,' whereas 'Magnesium,' 'Calcium' and 
'Osmolarity' were set to 'zero'. It was also assumed that in 
the immediate aftermath of engulfment by the macrophage, 
the SPI-1 secretion system of S. typhimurium and its main 
regulators (viz., HilD, HilC, RtsA and HilA) would still be 
in an active state. Thus, the state of these nodes were 
set to 'one' in the model corresponding to the intra- 
macrophageal environment. In summary, the initial condi- 
tions for the simulation were defined in order to reflect the 
in vivo state of human distal intestine and macrophage dur- 
ing the infection by S. typhimurium. 
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Additional file 3: Transcription factor binding sites predicted by 
Tfsitescan. Transcription factor binding sites predicted by 'Tfsitescan" in 
the upstream regions of Type VI Secretion System genes (sc/'S, vrgS and 
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pertaining to homology modeling of the proteins FlrC and YfhA and 
results of superposition of the two modelled structures. 

Additional file 6: Pfam Domains identified in FlrC and YfhA. 

Additional file 7: Homology modeling of three domains of FlrC and 
YfhA. Details pertaining to homology modeling of the 3 distinct 
domains from the proteins FlrC and YfhA, and the results obtained 
through pair-wise superposition of the corresponding domains from the 
two different proteins. 

Additional file 8: Stochastic state transitions. Analysis of state 
transitions obtained by simulating the Boolean model allowing 
asynchronous updates. 

Additional file 9: Boolean model simulation with mutated yfhA. 

Transition of states obtained through Boolean model simulation with 
mutated yfhA. 

Additional file 10: Logical rules for the Boolean model. Logical rules 
defined for activation of each of the 34 nodes in the Boolean model. 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

SSM and HR conceived the idea and design of the methodology. HR carried 
out initial analyses on a smaller Boolean network. CD, AD and SSM 
generated a comprehensive regulatory network based on the available 
literature. CD and AD carried out detailed analysis and simulation of the 
Boolean network. CD, AD, and SSM analysed all the results and drafted the 
manuscript. All authors read and approved the manuscript. 

Acknowledgements 

We thank Tungadri Bose, Mohammed Monzoorul Haque and Tarini Shankar 
Ghosh for helpful discussion and comments. 

Author details 

^io-Sciences R&D Division, TCS Innovation Labs, Tata Consultancy Services 
Ltd., 54-B, Hadapsar Industrial Estate, Pune 411013, Maharashtra, India. 
2 Present address: Novartis Healthcare Pvt. Ltd., #6 Raheja Mindspace, 
Hitec-city, Hyderabad 500081, India. 

Received: 16 July 2013 Accepted: 14 September 2013 
Published: 30 September 2013 

References 

1. Lucas RL, Lee CA: Unravelling the mysteries of virulence gene regulation 
in Salmonella typhimurium. Mol Microbiol 2000, 36:1024-1033. 

2. Marcus SL, Brumell JH, Pfeifer CG, Finlay BB: Salmonella pathogenicity 
islands: big virulence in small packages. Microbes Infect 2000, 2:145-156. 

3. Cirillo DM, Valdivia RH, Monack DM, Falkow S: Macrophage-dependent 
induction of the Salmonella pathogenicity island 2 type III secretion system 
and its role in intracellular survival. Mol Microbiol}998, 30:175-188. 

4. Hensel M: Salmonella pathogenicity island 2. Mol Microbiol 2000, 36:1 01 5-1 023. 

5. Bustamante VH, Martinez LC, Santana FJ, Knodler LA, Steele-Mortimer O, 
Puente JL: HilD-mediated transcriptional cross-talk between SPI-1 and 
SPI-2. Proc Natl Acad Sci USA 2008, 105:14591-14596. 

6. Martinez LC, Yakhnin H, Camacho Ml, Georgellis D, Babitzke P, Puente JL, 
Bustamante VH: Integration of a complex regulatory cascade involving 
the SirA/BarA and Csr global regulatory systems that controls expression 
of the Salmonella SPI-1 and SPI-2 virulence regulons through HilD. Mol 
Microbiol 2011,80:1637-1656. 

7. Leung KY, Siame BA, Snowball H, Mok Y-K: Type VI secretion regulation: 
crosstalk and intracellular communication. Curr Opin Microbiol 201 1, 14:9-15. 



Additional file 1: Components of the gene regulatory network. List 
of the genes (and RNAs) included in the network representing the cross-talk 
among regulatory elements of Type III and Type VI Secretion Systems in 
Salmonella typhimurium. The synonym codes correspond to the strain S. 
typhimurium LT2. 

Additional file 2: YfhA binding sites. YfhA binding sites found in the 
upstream region and in the ORF of res B. YfhA is known to recognize an 18 
base pair long motif [TGTCN(10)GACA]. The results given below indicate the 
position of the probable binding sites (with respect to the transcription start 
site), and the number of mismatches allowed while considering a hit. 



Das et al. Gut Pathogens 2013, 5:28 
http://www.gutpathogens.eom/content/5/1/28 



Page 12 of 12 



8. Bernard CS, Brunet YR, Gueguen E, Cascales E: Nooks and crannies in type 
VI secretion regulation. J Bacterid 2010, 192:3850-3860. 

9. Parsons DA, Heffron F: sciS, an icmF homolog in Salmonella enterica 
serovar Typhimurium, limits intracellular replication and decreases 
virulence. Infect Immun 2005, 73:4338-4345. 

10. Altier C, Suyemoto M, Lawhon SD: Regulation of Salmonella enterica serovar 
typhimurium invasion genes by csrA. Infect Immun 2000, 68:6790-6797. 

1 1 . Ellermeier CD, Ellermeier JR, Slauch JM: HilD, HilC and RtsA constitute a 
feed forward loop that controls expression of the SPI1 type three 
secretion system regulator hilA in Salmonella enterica serovar 
Typhimurium. Mol Microbiol 2005, 57:691-705. 

12. Ganesh AB, Rajasingh H, Mande SS: Mathematical modeling of regulation 
of type III secretion system in Salmonella enterica serovar Typhimurium 
by SirA. In Silico Biol (Gedrukt) 2009, 9:S57-S72. 

13. Fass E, Groisman EA: Control of Salmonella pathogenicity island-2 gene 
expression. Curr Opin Microbiol 2009, 12:199-204. 

1 4. Teixido L, Carrasco B, Alonso JC, Barbe J, Campoy S: Fur activates the expression 
of Salmonella enterica pathogenicity island 1 by directly interacting with the 
hilD operator in vivo and in vitro. PLoS One 201 1 , 6:e1 971 1 . 

1 5. Soncini FC, Groisman EA: Two-component regulatory systems can interact to 
process multiple environmental signals. J Bacteriol 1996, 178:6796-6801. 

16. Garcia-Calderon CB, Casadesus J, Ramos-Morales F: Regulation of igaA and 
the Res system by the MviA response regulator in Salmonella enterica. 

J Bacteriol 2009, 191:2743-2752. 

17. Majdalani N, Gottesman S: The Res phosphorelay: a complex signal 
transduction system. Annu Rev Microbiol 2005, 59:379-405. 

18. Wang Q, Zhao Y, McClelland M, Harshey RM: The RcsCDB signaling system 
and swarming motility in Salmonella enterica serovar typhimurium: dual 
regulation of flagellar and SPI-2 virulence genes. J Bacteriol 2007, 
189:8447-8457. 

19. Albert R: Boolean modeling of genetic regulatory networks. In Complex 
networks. Edited by Ben-Nairn E, Frauenfelder H, Toroczkai Z. Berlin 
Heidelberg: Springer; 2004:459-481. Lecture Notes in Physics, vol. 650. 

20. Bornholdt S: Boolean network models of cellular regulation: prospects 
and limitations. J R Soc Interface 2008, 5(Suppl 1):S85-S94. 

21 . Li F, Long T, Lu Y, Ouyang Q, Tang C: The yeast cell-cycle network is 
robustly designed. Proc Natl Acad Sci USA 2004, 101:4781-4786. 

22. Faure A, Naldi A, Chaouiya C, Thieffry D: Dynamical analysis of a generic 
Boolean model for the control of the mammalian cell cycle. 
Bioinformatics 2006, 22:e124-e131. 

23. Thakar J, Pilione M, Kirimanjeswara G, Harvill ET, Albert R: Modeling systems- 
level regulation of host immune responses. PLoS Comput Biol 2007, 3:e109. 

24. Giacomantonio CE, Goodhill GJ: A Boolean model of the gene regulatory 
network underlying Mammalian cortical area development. PLoS Comput 
Biol 2010, 6:e1 000936. 

25. Hegde SR, Rajasingh H, Das C, Mande SS, Mande SC: Understanding 
communication signals during mycobacterial latency through predicted 
genome-wide protein interactions and boolean modeling. PLoS ONE 
2012, 7:e33893. 

26. Garmendia J, Beuzon CR, Ruiz-Albert J, Holden DW: The roles of SsrA-SsrB and 
OmpR-EnvZ in the regulation of genes encoding the Salmonella typhimurium 
SPI-2 type III secretion system. Microbiology 2003, 149(Pt 9)2385-2396. 

27. Yamamoto K, Hirao K, Oshima T, Aiba H, Utsumi R, Ishihama A: Functional 
characterization in vitro of all two-component signal transduction 
systems from Escherichia coli. J Biol Chem 2005, 280:1448-1456. 

28. Reichenbach B, Gopel Y, Gorke B: Dual control by perfectly overlapping 
sigma 54- and sigma 70- promoters adjusts small RNA GlmY expression 
to different environmental signals. Mol Microbiol 2009, 74:1054-1070. 

29. Correa NE, Klose KE: Characterization of enhancer binding by the Vibrio 
cholerae flagellar regulatory protein FlrC. J Bacteriol 2005, 187:3158-3170. 

30. Syed KA, Beyhan S, Correa N, Queen J, Liu J, Peng F, Satchell KJF, Yildiz F, 
Klose KE: The vibrio cholerae flagellar regulatory hierarchy controls 
expression of virulence factors. J Bacteriol 2009, 191:6555-6570. 

31. Saini S, Slauch JM, Aldridge PD, Rao CV: Role of cross talk in regulating the 
dynamic expression of the flagellar Salmonella pathogenicity island 1 
and type 1 fimbrial genes. J Bacteriol 2010, 192:5767-5777. 

32. Levitt M: Accurate modeling of protein conformation by automatic 
segment matching. J Mol Biol 1992, 226:507-533. 

33. Finn RD, Tate J, Mistry J, Coggill PC, Sammut SJ, Hotz H-R, Ceric G, Forslund 
K, Eddy SR, Sonnhammer ELL, Bateman A: The Pfam protein families 
database. Nucleic Acids Res 2008, 36(Database issue):D281-D288. 



34. Feng X, Oropeza R, Kenney LJ: Dual regulation by phospho-OmpR of ssrA/ 
B gene expression in Salmonella pathogenicity island 2. Mol Microbiol 

2003, 48:1131-1143. 

35. Gonzalez AG, Naldi A, Sanchez L, Thieffry D, Chaouiya C: GINsim: a software 
suite for the qualitative modelling, simulation and analysis of regulatory 
networks. Biosystems 2006, 84:91-100. 

36. Wang M, Luo Z, Du H, Xu S, Ni B, Zhang H, Sheng X, Xu H, Huang X: 
Molecular characterization of a functional type VI secretion system in 
Salmonella enterica serovar Typhi. Curr Microbiol 201 1, 63:22-31. 

37. Bijlsma JJE, Groisman EA: The PhoP/PhoQ system controls the 
intramacrophage type three secretion system of Salmonella enterica. 
Mol Microbiol 2005, 57:85-96. 

38. Le D-H, Kwon Y-K A coherent feedforward loop design principle to sustain 
robustness of biological networks. Bioinformatics 201 3, 29:630-637. 

39. Mao F, Dam P, Chou J, Olman V, Xu Y: DOOR: a database for prokaryotic 
operons. Nucleic Acids Res 2009, 37(Database issue):D459-D463. 

40. Ghosh D: Object-oriented transcription factors database (ooTFD). Nucleic 
Acids Res 2000, 28:308-310. 

41 . Schwede T, Kopp J, Guex N, Peitsch MC: SWISS-MODEL: an automated 
protein homology-modeling server. Nucleic Acids Res 2003, 31:3381-3385. 

42. Zhu J, Weng Z: FAST: a novel protein structure alignment algorithm. 
Proteins 2005, 58:618-627. 

43. Ellermeier JR, Slauch JM: Fur regulates expression of the Salmonella 
pathogenicity island 1 type III secretion system through HilD. J Bacteriol 
2008, 190:476-486. 

44. Mizusaki H, Takaya A, Yamamoto T, Aizawa S: Signal pathway in salt- 
activated expression of the Salmonella pathogenicity island 1 type III 
secretion system in Salmonella enterica serovar Typhimurium. J Bacteriol 
2008, 190:4624-4631. 

45. Queiroz MH, Madrid C, Paytubi S, Balsalobre C, Juarez A: Integration host factor 
alleviates H-NS silencing of the Salmonella enterica serovar Typhimurium 
master regulator of SPI1, hilA. Microbiology 201 1, 157(Pt 9)2504-2514. 

46. Lim S, Yun J, Yoon H, Park C, Kim B, Jeon B, Kim D, Ryu S: MIc regulation of 
Salmonella pathogenicity island I gene expression via hilE repression. 
Nucleic Acids Res 2007, 35:1822-1832. 

47. Bowden SD, Rowley G, Hinton JCD, Thompson A: Glucose and glycolysis are 
required for the successful infection of macrophages and mice by 
Salmonella enterica serovar typhimurium. Infect Immun 2009, 77:31 17-3126. 

48. Miao EA, Miller SI: A conserved amino acid sequence directing 
intracellular type III secretion by Salmonella typhimurium. Proc Natl Acad 
Sci USA 2000, 97:7539-7544. 

49. Lee AK, Detweiler CS, Falkow S: OmpR regulates the two-component 
system SsrA-ssrB in Salmonella pathogenicity island 2. J Bacteriol 2000, 
182:771-781. 

50. Linehan SA, Rytkonen A, Yu X-J, Liu M, Holden DW: SlyA regulates function 
of Salmonella pathogenicity island 2 (SPI-2) and expression of SPI-2 
-associated genes. Infect Immun 2005, 73:4354-4362. 

51 . Mair SM, Nairz M, Bellmann-Weiler R, Muehlbacher T, Schroll A, Theurl I, 
Moser PL, Talasz H, Fang FC, Weiss G: Nifedipine affects the course of 
Salmonella enterica serovar Typhimurium infection by modulating 
macrophage iron homeostasis. J Infect Dis 201 1, 204:685-694. 

52. Ibarra JA, Knodler LA, Sturdevant DE, Virtaneva K, Carmody AB, Fischer ER, 
Porcella SF, Steele-Mortimer O: Induction of Salmonella pathogenicity 
island 1 under different growth conditions can affect Salmonella-host 
cell interactions in vitro. Microbiology 2010, 156(Pt 4):1 120-1 133. 

53. Eriksson S, Lucchini S, Thompson A, Rhen M, Hinton JCD: Unravelling the 
biology of macrophage infection by gene expression profiling of 
intracellular Salmonella enterica. Mol Microbiol 2003, 47:103-1 18. 

54. Craig M, Slauch JM: Phagocytic Superoxide Specifically Damages an 
Extracytoplasmic Target to Inhibit or Kill Salmonella. PLoS One 2009, 4:e4975. 

55. Gallois A, Klein JR, Allen LA, Jones BD, Nauseef WM: Salmonella 
pathogenicity island 2-encoded type III secretion system mediates 
exclusion of NADPH oxidase assembly from the phagosomal membrane. 
J Immunol 2001, 166:5741-5748. 

56. Albert I, Thakar J, Li S, Zhang R, Albert R: Boolean network simulations for 
life scientists. Source Code Biol Med 2008, 3:16. 



doi:1 0.1 1 86/1 757-4749-5-28 

Cite this article as: Das et al:. Understanding the sequential activation of 
Type III and Type VI Secretion Systems in Salmonella typhimurium using 
Boolean modeling. Gut Pathogens 2013 5:28. 



